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ABSTRACT 

The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, 
substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel 
PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This 
photometry indicates that a substantial mass of cold dust is present (> 3 x 10^ Mq) at temperatures significantly lower (20-28 K) than 
previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the 
cold gas traced by CO with a gas-to-dust ratio of 80-120. 
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. !^ ' 1 . Introduction that the dust continuum traces a substantial, cold component to 

S^ ^, ^ , ^ , . .... the ISM in these massive elliptical galaxies. However, the origin 

^ The cores of cluster of galaxies are very energetic regions with ^^ ^j^^ ^^^^ ^^^ ^^^ -^ -^ ^^-^^^^^ ^^.^ ^^. „ ^^^.j understood. 

rrt a high X-ray emissivity, particle density, cosmic ray flux, stellar ^, ,. . . ... 

■ ■ ■ density and AGN radiation. In this very hostile environment any . ^he limitations with the current observations of dust emis- 
dust grains are unUkely to survive for more than a few million ^1°" "i^^^e it difficult to establish an unambiguous dust mass as 
years due to the action of collisional sputtering (Dwek & Arendt *ey do not sample over the peak of the dust emission in the 
1992) unless they are shielded (Fabian et al. 1994). It is there- FIR^ The unprecedented sensitivity of Herschel (Pilbratt et al. 
fore somewhat surprising to find that dust continuum emission ^010) to FIR continuum offers the opportunity to accurately con- 
from the brightest cluster galaxies in the most rapidly cooUng !5^i" *^ f^" ^^^ ^P^'^^^"^ f'^f 'l^^t ^"^^^^1°" i" '^l^^'^^" ^^'^^'■ 
clusters being detected at sub-mm and MIR wavelengths (Edge ^he authors were awarded 140 hours of time in an Open Time 
et al. 1999, Egami et al. 2006, O'Dea et al. 2008). The presence ^^^ Project (PI Edge) to investigate the FIR line and continuum 
of cold molecular gas (Edge 2001 , Salome & Combes 2003) and Properties of a sample of 1 brightest cluster galaxies (BCGs) m 
dust absorption in HST imaging (McNamara et al 1996) implies well-studied cooling flow clusters selected on the basis of optical 
emission line and X-ray properties. The full goals of the project 



* Herschel is an ESA space observatory with science instruments are to observe at least five atomic cooling lines for each object 
provided by European-led Principal Investigator consortia and with im- that cover a range in density and temperature behaviour and ob- 
portant participation from NASA. tain a fully sampled FIR spectral energy distribution. In this pa- 
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per we present the Photodetector Array Camera & Spectrometer 
(PACS, Poglitsch et al. 2010) and Spectral and Photometric 
Imaging REceiver (SPIRE, Griffin et al. 2010) photometry for 
the three targets observed in the Science Demonstration Phase 
(SDP), Abell 1068 (z = 0.1386), Abell 2597 (z = 0.0821) and 
Zw3146 (z = 0.2906). In a parallel paper (Edge et al. 2010), we 
present the FIR spectroscopy for the first two of these clusters. 

The three clusters observed have quite contrasting multi- 
wavelength properties. Abell 1068 and Zw3 146 both have strong 
MIR emission (O'Dea et al. 2008, Egami et al. 2006) with a rela- 
tively bright CO detection (Edge 2001) and a weak central radio 
source (McNamara et al. 2004). A1068 lies just below the lumi- 
nosity threshold of a ULIRG (10'^ L©) and exhibits some contri- 
bution from an AGN (Crawford et al. 1999, O'Dea et al. 2008). 
On the other hand, Abell 2597 is a relatively weak MIR source 
(Donahue et al. 2007) with a weak CO detection (Salome, priv. 
comm.) and a powerful central radio source (Sarazin et al. 1995). 
The implied FIR luminosity of A2597 is a factor of around 30 
below that of A 1068 and, in addition, the fractional contribution 
from an AGN in the MIR is also lower in A2597. 



2. Observations 

We performed photometric imaging of A 1068, A2597 and 
Zw3146 with PACS and SPIRE. The data were reduced with 
the Herschel Interactive Processing Environment (HIPE) soft- 
ware version 2.3.1436 (Ott 2010). We used for both PACS and 
SPIRE the official scripts as presented by the PACS and SPIRE 
ICC teams during the Herschel SDP data processing workshop 
in December 2009. 



2.1. PACS Data 

The PACS photometric observations were taken in 
LargeScanMapping mode in all three bands of the pho- 
tometer, BS (70 jum), BL (100 yum) and R (160 fim) using the 
medium scan speed (20"s"'). The scan maps comprised 18 scan 
line legs of 4' length and cross-scan step of 15". Each obser- 
vation had a "scan" and an orthogonal "cross-scan" direction 
and we calibrated the corresponding data separately before 
combining them into a single map of 9 X 9 . The resulting maps 
have a resolution of 5.2", 7.7" and 12" at 70, 100 and 160 fim, 
respectively and are presented in the electronic version of this 
paper The PACS photometer performs dual-band imaging such 
that the BS and BL bands each have simultaneous observations 
in the R band so we have two sets of scans in the R band. 

We adopted the PACS Data Reduction Guideline to process 
the raw level-0 data to calibrated level 2 products and used the 
official script for PACS ScanMapping mode but with particular 
attention to the high pass filtering to remove "l/-\/7" noise. We 
choose to use the HighPassFilter method with a filter of 20 read- 
outs which will remove structure on all scales above 82". The 
target BCG and other bright sources in the field were masked 
prior to applying the filter. The size of the mask was chosen to 
be less than the filter size so as to minimize any left-over low- 
frequency artefacts under the masks. We used masks with a ra- 
dius of 15" for our sources. We tried varying the size for the 
filter from 10 to 30 readouts and the mask radius from 10-30" 
and found our results to not change significantly for these ranges 
in values. Finally the task 'photProject', was used to project the 
calibrated data onto a map on the sky in units of Jy pixel'. The 
"scan" and "cross-scan" maps were then averaged to produce the 
final coadded map. The PACS and SPIRE images are included 



in the electronic version of the paper. The spatial flux distribu- 
tion and flux densities of our target sources were investigated 
using cumulative flux curves. The spatial flux distribution for 
each of our three sources is consistent with that expected from a 
point source. Flux densities in the BS, BL and R band were ex- 
tracted using a 33" by 33" aperture centered on the BCG. Small 
aperture corrections were applied as outlined in the PACS Scan 
Map release note (PICC-ME-TN-035). Care was taken to cali- 
brate these derived flux densities to account for the known flux 
overestimation in the used HIPE version by factors 1.05, 1.09 
and 1.29 in BS, BL and R bands respectively. The absolute flux 
accuracy is within 10 % for BS and BL, and better than 20 % 
for R. These uncertainties are not believed to be correlated due 
to the BS and BL bands being taken at different times and the R 
band using a diff'erent detector 

2.2. SPIRE Data 

The SPIRE photometry was performed in the LargeScanMap 
mode with cross-linked scans in two orthogonal scan directions. 
The photometer has a field of view of 4 x 8 , which is ob- 
served simultaneously in three spectral bands, PSW (250 jum), 
PMW (350 nm) and PLW (500 //m) with a resolution of about 
18", 25" and 36", respectively. The resulting maps measure 
12 X 12 in size and are presented in the electronic version of 
this paper. 

We used the standard HIPE pipeline for the LargeScanMap 
observing mode and the naive map-maker. The pre-processed 
raw telemetry data were first subject to engineering conver- 
sion wherein the raw timeline data were converted to mean- 
ingful units, the SPIRE pointing product was created, deglitch- 
ing and temperature drift correction were performed, and maps 
were created, the units of which were Jy beam"'. Our targets 
are unresolved at the spatial resolution of SPIRE. We derived 
their flux densities by fitting the sources with the SPIRE point 
source response function. Care was taken to de-blend our tar- 
get from other nearby sources at the longer wavelengths, where 
the sources are most likely to be background to the cluster We 
account for the known flux calibration offset in the used ver- 
sion of HIPE by applying the following multiplicative calibra- 
tion factors 1.02, 1.05 and 0.94 to the derived flux densities in 
the PSW, PMW and PLW bands respectively (see Griffin et al. 
2010, Swinyard et al. 2010). We also performed aperture pho- 
tometry using the HIPE point-source extraction (PSE) tool but 
this method gives accurate results only for isolated point sources. 
At 350 jum and 500 //m, the BCGs in A2597 and Zw3146 are 
close to the detection limit and at the confusion limit of SPIRE 
making the PSE method of determining the fluxes unsuccessful. 
A 1068 has a relatively strong compact BCG in far infrared and 
so we performed the PSE to find that the flux estimates using 
AIPS and HIPE agree with each other to better than 5%. 

3. Results 

In the PACS photometry, A1068, A2597 and Zw3 146 have been 
detected in all three bands. For A1068, 70 and lOO/zm values are 
slightly less than the IRAS 60 and 100 //m measurements. This 
could be due to nearby sources that cannot be separated from the 
BCG in the much lower resolution IRAS observations but no suf- 
ficiently bright source is visible in our PACS imaging. There is a 
large difference between the Spitzer MIPS 70//m flux (Quillen et 
al. 2008) and our PACS 70yum flux, the PACS flux being a factor 
1.7 lower than the MIPS flux. In the case of Zw3146 the MIPS 
and PACS 70 /im fluxes also differ with the PACS value being 
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Table 1. Log of Herschel Observations. 



Cluster 



A 1068 



0.1386 



A2597 



0.0821 



Zw3146 0.2906 



Instrument 


A 


Obsid 


Flux 




(jj.m) 




(mJy) 


PACS 


70 


1342187051 


542+6 


PACS 


100 


1342187053 


757+6 


PACS 


160 




769±4 


SPIRE 


250 


1342187321 


376±6 


SPIRE 


350 




135±6 


SPIRE 


500 




56±8 


SCUBA 


450 




39±13 


SCUBA 


850 




5.3±1.1 


Spitzer 


24 




74.5+2.0 


Spitzer 


70 




941+30 


IRAS 


60 




577+52 


IRAS 


100 




958+144 


PACS 


70 


1342187118 


57+5 


PACS 


100 


1342187120 


67+7 


PACS 


160 




86+4 


SPIRE 


250 


1342187329 


30+6 


SPIRE 


350 




15+6 


SPIRE 


500 




< 16 


SCUBA 


850 




14.5+2.3 


Spitzer 


24 




2.1+0.2 


Spitzer 


70 




49+6 


Spitzer 


160 




52+3 


PACS 


70 


1342187043 


94+6 


PACS 


100 


1342187045 


150+6 


PACS 


160 




139+5 


SPIRE 


250 


1342187326 


81+6 


SPIRE 


350 




30+6 


SPIRE 


500 




< 16 


SCUBA 


450 




<48 


SCUBA 


850 




6.6+2.6 


Spitzer 


24 




4.1+0.4 


Spitzer 


70 




68+14 


Spitzer 


160 




157+35 



Notes. The Spitzer data are from Quillen et al. (2008), Donahue et al. 
(2007, priv. comm.) and Egami et al. (2006). The SCUBA data are from 
Edge (priv. comm.), Zemcov et al. (2007) and Chapman et al. (2002). 



a factor 1.4 larger than the MIPS value (Egami et al. 2006). For 
A2597 the PACS fluxes diff^er from the Spitzer 70 and 160 fim 
fluxes reported by Donahue et al. (2007). Part of this difference 
was resolved when the MIPS 70 j^m data were re-analysed and 
found to be a factor of two too high (Donahue, priv. comm.). 
The differences observed between the PACS and Spitzer fluxes 
require further investigation. In the SPIRE photometry, A 1068 
is detected in afl three SPIRE bands. A2597 and Zw3146, while 
clearly detected in PSW and PMW bands, have a 1-2 cr detec- 
tion in the PLW band. Table 1 gives the photometric results for 
the three galaxies, with 2cr upper-limit for A2597 and Zw3 146 in 
PLW. Figure 1 presents the radio to optical spectral energy dis- 
tributions (SEDs) for the three targets. These plots show the sig- 
nificant variation in the relative radio-FIR-optical contributions 
for each of our galaxies. Here we focus on the sub-mm/MIR dust 
emission as sampled by PACS and SPIRE photometry, comple- 
mented by published Spitzer and IRAS measurements. 

We fit the SEDs of the dust emission using black bodies 
modified with a dust emissivity index, /3. The FIR-MIR slopes 
of our sources require the presence of at least two dust com- 
ponents. Previous studies of star-forming galaxies have indeed 
established that a single modified black body (MBB) is inade- 
quate to account for the observed dust emission (Wiklind 2003). 
Hence, our model for the SEDs consists of two MBBs with the 



dust emissivity index for each fixed to /3=2 and a mass absorp- 
tion coeflicient, /Cdv, of 2.5 m^ kg~' at lOOyum. 

For A 1068 we fit the 24-850yum emission. For A2597 and 
Zw3146 the SCUBA 850/im detections have been removed and 
we fit only the 24-350/im range. In the case of A2597, this is 
due to the unknown amount of radio contamination at 850/im. In 
the case of Zw3146 the BCG is blended with strong background 
source at 850/im (Chapman et al. 2002). The data are weighted 
in the fit inversely to the square of their error. The resulting fits 
are shown in Figure 1 . The derived dust temperatures and total 
FIR luminosities for each source are listed in Table 2. 

The results in Table 2 indicate that at least two dust compo- 
nents, one at 20-25 K and one at 50-60 K, are present in all three 
sources. The FIR emission is much stronger relative to the opti- 
cal in A1068 and Zw3146 as compared to A2597. The SEDs of 
A1068 and Zw3 146 resemble those of strongly star-forming sys- 
tems and, based on the total FIR luminosity derived here, we find 
star formation rates (SFR)of 6Oand44M0yr"' in these two sys- 
tems using the Kennicutt (1998) conversion factor. For A2597 a 
much more modest SFR of 2 M© yr"' is inferred. These values 
are comparable to SFRs derived from Ha line and/or UV contin- 
uum emission given the uncertainties of these tracers. However, 
the SFR values derived from Spitzer data are higher for A 1068 
and Zw3146. The difference for A1068 is the most pronounced 
and can be directly attributed to the stronger AGN contribution 
in this object (Quillen et al. 2008) which boosts the 24/i flux 
compared other comparable sources. Therefore, when the to- 
tal FIR luminosity is derived from the 15/im flux infered from 
Spitzer it will be overestimated. The value for Zw3146 from 
Egami et al. (2006) is higher than ours as their fit includes the 
SCUBA 850/im point from Chapman et al. (2002) which appears 
to be overestimated on the basis of our SPIRE data. 

The gas to dust ratio is found to be between 80 and 140 (see 
Table 2). Gas temperatures can be inferred from CO measure- 
ments (Edge 2001, Salome & Combes 2003). These estimates 
infer gas temperatures of 25-40 K thus implying that the gas 
and dust share a common environment and are potentially co- 
located in the denser regions of cold, molecular gas clouds. We 
have attempted to determine how much extended emission is 
present from our highest spatial resolution PACS lOfim image 
but we find no evidence for more than 10% additional flux be- 
yond a point source. Clearly these limits will improve with a 
better characterisation of the instrument but we believe that we 
can conclude that the dust emission in our targets has an extent 
comparable to that the bulk of the CO emitting gas and optical 
emission lines (< 5" or 5-20 kpc). 

4. Discussion and conclusions 

Our initial Herschel results confirm the presence of the striking 
dust emission peak expected from the observations at sub-mm 
(Edge et al. 1999, Chapman et al. 2002) and MIR (Egami et al. 
2006, O'Deaetal. 2008). 

The star formation rates derived from the full-sampled FIR 
SED are comparable to those derived from Spitzer 24//m fluxes 
apart from A1068, which has the strongest contribution from 
an AGN so hot dust dominates to the 24//m flux. However, in 
the sub-mm the contribution from the radio continuum from an 
active nucleus must be correctly accounted for before any dust 
mass can be estimated from the 850/im flux. In the case of A2597 
here and A2390 in Edge et al. (1999), the presence of a powerful 
radio source appears to contribute to the SCUBA 850/im flux. 

While it is difficult to draw any general conclusions from just 
three BCGs, we note with interest that the ratio of dust mass to 
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Table 2. Summary of results and other cluster properties. 



Cluster 



A 1068 



A2597 



Dust Temperatures 


24±4K 


21±6K 


23±5K 




5Tl^K 


48+^'K 


53!fK 


Cold Dust Mass 


5.1x10* Mo 


2.3x10' Mo 


5.4xlO'*Mo 


Warm Dust Mass 


3.9x10* Mo 


2.9x10' Mo 


1.9x10" Mo 


Total FIR Luminosity 


3.5x10" Lo 


8.8x10' Lo 


2.5x10" Lo 


Star Formation Rate 


60+20 Mo yr-' 


2+lMoyr-' 


44+14 Mo yr-' 


SFR Spitzer 


188 Mo yr-' 


4 Mo yr-' 


70±14Moyr-' 


SFR optical/UV 


20-70 Mo yr-' 


10-15 Mo yr-' 


47±5 Mo yr-' 


CO gas mass 


4.1x10'" Mo 


2.0x10' Mo 


7.7x10'" Mo 


Ho- Slit Luminosity 


8x10^" ergs-' 


3x10^" ergs-' 


3xlO''2ergs-' 



Sarazin, C. L., Burns, J. O., Roettiger, K., et al. 1995, ApJ, 447, 559 
Swinyard, B., Griffin, M., Ade, P., et al. 2010, A&A, this volume 
====mm'G. B., O'Dea, C. P., Peck, A. B., et al. 1999, ApJ 512, L27 
WiHind, T. 2003, ApJ, 588, 736 
__Zemcov, M., Borys, C, Halpem, M., et al. 2007, MNRAS, 376, 1073 



Notes. The Spitzer SFR values are from O'Dea et al. (2008), Donahue 
et al. (2007) and Egami et al. (2006). The Optical/UV SFR values are 
from McNamara et al. (2004), Donahue et al. (2007) and Egami et al. 
(2006). The CO gas masses are from Edge (2001) and Salome (priv. 
comm.) and the Ho- slit luminosities are from Crawford et al. (1999). 



CO-derived gas mass is constistent for all three within a factor 
of five. If the dust were mostly generated through dust ejection 
from evolved stars then the dust mass should closely correlate 
with the total stellar mass. However, our three galaxies have very 
similar optical/NIR absolute magnitudes. So, unless the ejected 
dust were "captured" by the cold gas clouds protecting it from 
X-ray sputtering, this suggests that the apparent correlation be- 
tween the molecular gas and dust masses arises from a direct 
connection between the gas reservoir and star formation. 

These results are a limited example of those to come in the 
very near future from Herschel as there are two other Open Time 
Key Projects (PI Egami and Smith) that are targetting a total of 
70 clusters that cover a broad range of BCG properties so the 
wider context of these initial observations can be determined. In 
particular, the amount of dust present in more quiescent BCGs 
and other massive cluster ellipticals will be important in assess- 
ing how much of the dust seen in cool core BCGs orginates from 
the underlying stellar population. 
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Fig. 1. Spectral energy distributions for A1068 (top), 
A2597 (middle) and Zw3146 (bottom) including Herschel 
PACS/SPIRE (blue/red symbols), Spitzer, Radio, NIR photom- 
etry from 2MASS and optical photometry from SDSS. To ac- 
count for absolute flux uncertainties we have set the following 
errors on the fluxes derived from the various instruments (unless 
the quoted error is larger than this); PACS BS/BL 10%, PACS 
R 20%, SPIRE 15%, Spitzer 10%, SCUBA 450 jjm 30% and 
SCUBA 850 yum 20%. The model fit to the sub-mm/FIR/MIR 
data is shown by the black solid line. Only filled symbols have 
been used in the fit. The two modified blackbodies making up the 
model are shown by the black long dash and dash-dot lines. For 
A2597 we also show two VLBI measurements (black crosses) 
of the BCG core at 1.3 and 5 GHz (Taylor et al. 1999). These 
points show that the BCG has a strong, inverted radio core. 
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A1068 



A2597 



Zw3146 



Fig. 2. Colour images from the three PACS bands (BS, BL and R in the blue, green and red channels) for the three clusters within 
radius of 2.5' of the BCG The top row are images combined in their original resolution and the bottom row are the images combined 
with a common smoothing of 12" to match resolution. 
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A1068 



A2597 



Zw3146 



Fig. 3. Colour images from the three SPIRE bands (PS W, PMW and PLW in the blue, green and red channels) for full field covered 
for the three clusters covering approximately 12' x 12'. The top row are images combined in their original resolution and the bottom 
row are the images combined with a common smoothing of 36" to match resolution and clipped to remove areas of low exposure. 
The BCG is at the centre of the image and in A2597 and Zw3146 is the bluest object present (see text). 



